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After a brief introduction to the Heavy-Quark Effective Theory (HQET), I review the extraction of the Isgur- 
Wise function from lattice QCD calculations of the matrix elements for semi-leptonic decays of heavy-light pseu- 
doscalar mesons both into pseudoscalar and into vector mesons. This work is beginning to test the heavy-quark 
spin-flavour symmetries around the charm mass and to indicate the size of 0(l/m c ) corrections. An alternative 
approach to put the HQET on the lattice offers the prospect of computing the Isgur-Wise function directly. 



1. INTRODUCTION 

1.1. HQET 

The motivation behind Isgur-Wise function cal- 
culations is to be able to extract CKM matrix ele- 
ments from experimental measurements of heavy- 
quark decays in a model-independent way. Of 
particular interest are the semi-leptonic decays 
B — > Dlv and B — > D* lv, from which it is possi- 
ble to extract V c j. 

I shall consider only heavy-light mesons in 
which the heavy-quark mass, tuq , is much greater 
than Aqcd- Then the four- velocity of the heavy 
quark is essentially the same as the four-velocity 
of the meson and, in its rest frame, the heavy 
quark appears to the light degrees of freedom as 
a static colour source. The spin interaction be- 
tween the heavy quark and the light degrees of 
freedom is suppressed by the heavy-quark mass. 
So, if a weak current changes the heavy quark 
into a heavy quark of a different flavour (mass), or 
flips its spin, without changing its velocity, then 
the light degrees of freedom see no change. This is 
the reason for the spin-flavour symmetries of the 
HQET. If the weak current changes the velocity 
of the heavy quark, then the light degrees of free- 
dom have to adjust themselves, and this results in 
a form-factor suppression. The overlap between 
the state of the light degrees of freedom in the 
presence of a colour source moving with veloc- 
ity v, and that in the presence of a colour source 
moving with velocity v' , is the Isgur-Wise func- 
tion [1]. This is a fundamental quantity in QCD, 



which can only be determined non-perturbatively. 

Since the velocity of the heavy quark is con- 
served by soft QCD processes, it may be used to 
project out the upper spinor components of the 
heavy-quark field, Q, 



1 + y 



Q = e -im Q vx Ht 



i>H = H. 



(1) 



The HQET is obtained by integrating out the 
lower components [2]. There is a different field 
H for each velocity [3]. A semi-leptonic decay 
involves a transition between a heavy quark with 
velocity v and another with velocity v' so, at lead- 
ing order in I/tuq, the Lagrangian is 



HQET 



Hiv ■ DH + H'iv' ■ DH'. 



(2) 



1.2. Semi-leptonic form factors 

The matrix elements of the vector and axial 
vector currents between pseudoscalar (P) and 
vector (P*) meson states may be expressed in 
terms of six form factors, hi = hi(u>): 

(P' t v'\Vn\P, v) = (v + v')^h+ + (v- v')^h_ (3) 



{P*, v',e\V^\P, v) = ie^ pa e 



(4) 



(P*,v',e\A ll \P,v) = (w 



v^h A2 +v'^h Ai \ (5) 



where v and v' are the velocities of the incoming 
and outgoing mesons, u> = v ■ v' , and I have used 
a mass-independent normalisation for the states. 
In the HQET, the spin-flavour symmetry imposes 
constraints on the form factors. Eq. (1) gives 
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(v - v')^H'^H = 



h_(uj) = Q. (6) 



The time component of the vector current gener- 
ates the flavour symmetry: 



P',v\ / d 3 xH'y°H\P,v 



(P',v\P',v) 
h+(l) = l. 



(7) 



Flipping the spin of the heavy quark changes a 
pseudoscalar meson into a vector meson. The full 
consequences of the heavy-quark symmetry are [4] 

(M',v'\H'TH\M, v) = -£(w)Ti[M'TM] (8) 

l + j / -75 for \P,v) 



M 



£ ioi\P*,v,£) 



(9) 



so that, in the limit tuq — > oo, 

h+ = h v = h Al = h Ai = t{w), = i, (io) 



0. 



(11) 



is the Isgur-Wise function. It is normalised 
to one at zero velocity transfer (zero recoil). 

1.3. Renormalised Isgur-Wise function 

The above results are only exact in the limit of 
exact spin-flavour symmetry. We want the corre- 
sponding results in QCD, when the heavy-quark 
mass is large but finite. QCD and the HQET con- 
tain the same long-distance physics, but differ at 
short distances, where QCD resolves the dynam- 
ics of the heavy quark through the high momen- 
tum that can flow in internal loops. These short- 
distance corrections are included by expressing a 
current in QCD as an expansion in terms of local 
operators in the HQET: 

J^ CD ~ £ Q(M)^(M) HQET + 0(l/m Q ). (12) 

i 

Although in QCD the vector and axial-vector 
currents are not renormalised, the correspond- 
ing HQET currents are not symmetry currents 
and so depend on the renormalisation scale /j,. In 
Eq. (12) this dependence is cancelled by that of 
the Wilson coefficients, C;. 

The Wilson coefficients have been calculated to 
next-to-leading order by Neubert [5]. This calcu- 
lation is complicated by the fact that there are 



two heavy-quark scales: raj and m c . Neubert's 
review [6] provides a complete discussion. He 
shows that the Wilson coefficients factorise: 



d(u),fj) = Ci(m b ,m c ,uj)K(uj, fj) 



(13) 



where K contains all the dependence on n, is in- 
dependent of the spin structure, and 



K{l,n) = l. 



(14) 



Therefore, up to corrections of 0{l /rag), the ma- 
trix element of a heavy-quark current in QCD is 
given in terms of HQET matrix elements by 

(M',v'\J\M,v) = J2 C i( M ' ,v'\H'TiH\M,v) 

i 

= -6en^QTr[A^'r^], (15) 



where 

£renH = £(w , fJ.)K (w , fj.) 



(16) 



is the renormalisation-group invariant Isgur-Wise 
function, which is also normalised to one at zero 
recoil, because of Eq. (14). 

The short-distance corrections break the 
heavy-quark symmetry so that, from Eq. (15), the 
naive relations between the form factors, Eq. (10) 
and (11), are replaced by 



'-(c 2 + c 3 ) 



'-(c 2 -c 3 )C, 



Cl £ r en 
V5 



= ^2 £ren 

h A3 = (C\ + cl) & 



(17) 

(18) 

(19) 
(20) 
(21) 
(22) 



Neubert [5,6] gives expressions for Ci and Cf as 
functions of u>, raj and m c . He finds that, for the 
physical b and c masses, the next-to-leading-order 
corrections can be as large as 15%. This result is 
important, because we extract estimates of £ ren 
from measurements of the form factors in QCD 
using Eq. (17)-(22), i.e., 
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Cren ~ -y — ^rad corr- 



(23) 



It is useful to parametrise the Isgur-Wise func- 
tion by means of its slope at zero recoil, £'(1) = 
— p 2 . Most models agree in the kinematically ac- 
cessible region, 1 < u> < 1.6, and so I will define 
p 2 using the Bauer-Stech-Wirbel (BSW) model: 



■ exp 



"(2P 2 - 1) 



(24) 



QCD sum rules give bounds on p 2 [6]: 

Bjorken: p 2 > 0.25 (25) 



Voloshin: p 2 < 0.25- 



m M - mQ 

-^min 



(26) 



1.4. Power corrections 

At order I/tuq there are new universal func- 
tions (see [6]). Lattice calculations are not yet 
precise enough to disentangle these. The most 
useful statement is Luke's Theorem [7]: 0(1 /itiq) 
corrections to meson-decay matrix elements van- 
ish at zero recoil. This only protects 



h + (l) = ^C i (l) + 0(l/m 2 Q ) 



h Al (l) = Ct(l) + 0(l/m 2 Q ) 



(27) 



(28) 



because the other form factors are multiplied by 
terms which vanish anyway at u> = 1. Using QCD 
sum rules for b — > c, Neubert [6] estimates the 
power corrections to be around 3% for h + and 
h,A 1 , throughout the kinematic range. In con- 
trast, the correction to the relation between hy 
and £ is around 30%. A challenge for lattice cal- 
culations, with important phenomenological im- 
plications, is the determination of ratios such as 
h-v/h-A! and h + /hA 1 - 

2. RESULTS FROM LATTICE QCD 

Bernard et al. [8,9] were the first to report 
results for the Isgur-Wise function extracted 
from form factors computed in quenched lat- 
tice QCD. UKQCD has now also reported re- 
sults [10]. Both groups have presented further 
analyses at this conference [11-15]. The param- 
eter values used are given in Table 1. Both 



Table 1 

QCD simulations using 16 3 and 24 3 lattices 
(3 q-^GeV) ro P (GeV) ro g - g (GeV) 
Bernard et al: Wilson action, ~ 20 configs 



2.0 
3.2 



1.6-2.3 
1.6-2.7 



0.6-1.0 
0.6-1.0 



UKQCD: SW action, 36 and 60 configs 



6.0 
6.2 



2.0 
2.7 



1.6-2.1 
1.5-2.4 



0.5-0.8 
0.5-0.8 



groups use heavy-quark masses around charm and 
light-quark masses around strange. UKQCD uses 
the 0(a)-improved Sheikholeslami-Wohlert (SW) 
fermion action [16], whereas Bernard et al. use 
the Wilson fermion action. 

2.1. Pseudoscalar form factors: P — > P' ' Iv 

For elastic scattering {juq = rriQi), in QCD, 
because of current conservation, 



L elastic 



M = o, 



Lelastic 



(l) 



(29) 



to all orders in 1/rriQ. It follows that, close to 
zero recoil, the 0(1 /itiq) differences between h + 
and £ should be small, and that ^racTcorr should 
provide a good determination of £ ren - Assum- 
ing that heavy-quark symmetry applies to b — > c 
transitions, this opens the possibility of a model- 
independent determination [6] of V c j from exper- 
imental measurement of 



dT(B -> D*lv) 



known factorx£ 2 (w)|V ci) | 2 .(30) 



The matrix elements of the lattice vector cur- 
rent between pseudoscalar states of various mo- 
menta are obtained in the usual way from 3- 
point correlation functions, using amplitudes and 
energies from fits to 2-point functions. To fix 
the normalisation of the lattice vector current 
(ZyV^ = V M ), both groups divide by the for- 
ward matrix element of Vo at zero momentum, 
which in the continuum is 2. This cancels Zy 
and any field normalisation, e.g., as is required 
for Wilson fermions at large quark mass [9] . This 
approach should tend to reduce the effect of lat- 
tice artefacts which render the normalisations un- 
certain. UKQCD data for the SW action show 
that Zy is insensitive to the light-quark mass, or 
to whether momentum or w/12 is used, but it 
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0.06 — 



V.V 

Figure 1. Isgur-Wise function obtained by 
Bernard et al. from h + (crosses) compared with 
ARGUS data (circles). 

does depend on the heavy-quark mass. Thus, the 
procedure is unambiguous for elastic matrix ele- 
ments, given that h + is normalised at zero recoil. 

The estimate of the Isgur-Wise function ob- 
tained by Bernard et al. [9,11] is shown in Fig- 
ure 1, along with experimental data from AR- 
GUS [17]. The lattice results are from a mixture 
of (3 and tuq values, at fixed light-quark mass 
around m s . No radiative corrections are applied. 

The corresponding results from UKQCD [10, 
12] are shown in Figure 2. Here the data are 
at fixed (3 (6.2), fixed heavy-quark mass (charm), 
have been chirally extrapolated, and have had ra- 
diative corrections applied as in Eq. (23). A fit of 
Eq. (24) to the lattice data gives 



1.4±2±4 
1.2 



+ 7 
- 3 



Bernard et al. 
UKQCD 



(31) 



where some estimate of systematic errors has 
been included by both groups. The corresponding 
BSW functions are shown as solid curves. Using 
these to extrapolate the experimental data to zero 
recoil gives 



\V C 



cb\ 



0.044±7±5 
0.038 



+ 2+8 
-2-3 



Bernard et al. 
UKQCD 



(32) 



where the first error is experimental and the sec- 
ond is from the lattice (the B lifetime is taken 



0.02 — 




Figure 2. Isgur-Wise function obtained by 
UKQCD from h + , compared with ARGUS data. 

to be 1.49ps). Systematic errors apart, the lat- 
tice calculations are quite precise, in comparison 
with the experimental data, and the agreement 
between them is good. 

I will now turn to the systematics. The form 
factor h + is obtained from the elastic matrix ele- 
ment by setting /i_ = (Eq. (29)). This is true in 
the continuum, but not necessarily so on the lat- 
tice. UKQCD has tested this [12] and finds that, 
at (3 = 6.2, the measured values of Zi_ are con- 
sistent with zero to within 2a, provided that the 
maximum momentum in any external particle is 
less than or equal to \/Z{tt/12). 

Also reported [12] are the results of a more re- 
fined analysis of the (3 = 6.2 UKQCD data, using 
a correlated fit to all non-zero matrix elements 
and more momentum combinations. This sug- 
gests that p 2 decreases slightly with light-quark 
mass to a value of O.9I3 in the chiral limit, al- 
though the data are consistent with a constant 
value around 1.2. Bernard et al. [9] also note this 
tendency, but their data are rather noisy. 

At (3 = 6.0, UKQCD determines h + and Zi_ 
at fixed u> > 1 and at fixed mp, as an expansion 
in the hopping parameter for the b quark [13]. 
Although the slight dependence of Zy on tuq 
is ignored in this preliminary analysis, UKQCD 
finds that Zi+ ra d corr is approximately indepen- 
dent of ms, so that power corrections are small, 
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l.o 1.1 1.2 1.3 

v.v' 

Figure 3. Isgur-Wise function, normalised to 
Z^ 1 , obtained from at fixed heavy-quark 

mass (circles), and from different heavy-quark 
masses at maximum recoil (crosses). 

as expected, indicating that the heavy-quark 
flavour symmetry holds at around the charm mass 
(within large statistical errors), whereas, h_/h + 
shows a stronger dependence on mj than the size 
of the short-distance and power corrections esti- 
mated by Neubert [6] would indicate. If system- 
atic errors can be controlled, such measurements 
will be important predictions for the lattice. 

2.2. Vector and axial-vector form factors: 

P -H- P*ll> 

So far, only UKQCD has reported results for 
the Isgur-Wise function extracted from matrix el- 
ements for the decay of a pseudoscalar into a vec- 
tor meson [14,15]. This preliminary analysis is 
providing the first tests of the heavy-quark spin 
symmetry in semi-leptonic decays. 

Although the lattice vector current may be nor- 
malised using (P,v\Vfj,\P,v) = 2u M , the normali- 
sation of the axial-vector current, ZaA 1 ^ = A^, 
must be determined separately. The UKQCD 
analysis assumes that Za is an overall constant, 
and shows that this assumption is supported 
(within relatively large statistical errors) by esti- 
mating Z^ 1 = Ha 1 (l)rad corr [14]. The further un- 
certainty in the normalisation of the quark fields 
is much reduced by using the SW action [18]. 



UKQCD preliminary 

k q = k v = 0.129 - 
K q = 0.14144 

: * * : 


'- X ZA(hv/h Al ) r ad corr 




z — * — # 

- ° ZA(h + A Al )rad corr C 


e — 



1.0 1.1 1.2 1.3 

v.v' 

Figure 4. Ratios of estimates of the Isgur-Wise 
function. 

Figure 3 [15] shows the Isgur-Wise function ex- 
tracted from the axial-vector form factor. Data 
are included at fixed equal heavy-quark masses, 
rriQ = rriQi , and various momenta, and at fixed 
momentum transfer (g 2 = 0) and various heavy- 
quark masses, tuq ^ rriQi , for which 

I 2 rn 1 ( mp , mp * \ f**\ 

"{V =°)=o + • ( 33 ) 

2 \tnp. mp J 

The agreement between the two sets of data pro- 
vides further evidence for heavy-quark flavour 
symmetry at around the charm mass. 

Evidence for spin symmetry at the charm mass 
is given by the fact that the ratio h + /hA 1 in Fig- 
ure 4 is approximately constant and close to 1. 
Both form factors in this ratio are protected by 
Luke's Theorem. In contrast, the ratio of the vec- 
tor to axial- vector form factors, also shown in Fig- 
ure 4, is significantly different from 1, although 
roughly constant [15]. This supports the QCD 
sum rule result that hy suffers large 0(1 /itiq) 
corrections. These appear to be roughly propor- 
tional to £ and around 40% at the charm quark 
mass. The size of these corrections and their 
dependence on the heavy-quark masses [15] are 
roughly consistent with the QCD sum rule esti- 
mate of a 30% effect for b — > c [6]. 
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3. LATTICE HQET 

The Isgur-Wise function may also be calculated 
by formulating the HQET, Eq. (2), directly on a 
lattice. Mandula and Ogilvie [19] solve a lattice 
version of the heavy-quark propagator equation, 

iv ■ DS(x,y;v) = S(x,y), (34) 

numerically, by forward recursion, whereas Agli- 
etti [20] expands S as a power series in the veloc- 
ity v, obtaining £'(1) at 0(v 2 ). 

At this conference, Mandula and Ogilvie pre- 
sented results for obtained by combining 
S with Wilson light-quark propagators at (3 = 
5.7 [21]. They encounter problems with noise 
and find it essential to smear the heavy-quarks 
in order to see a signal even at short time separa- 
tions. Consequently, they are unable to explore 
whether their data is asymptotic. They have the 
advantage that the forward recursion is numeri- 
cally easy and that v is a free parameter, so that 
many values of u> can be sampled. The Isgur- 
Wise function they obtain is consistent with the 
estimates from lattice QCD and with experiment, 
within quite large errors, although the renormali- 
sation factor has not yet been included [20]. From 
a quadratic fit they obtain £'(1) = —0.95. 

It will be important for understanding the size 
of power corrections to obtain a more precise es- 
timate by this method. 

4. CONCLUSIONS 

The Isgur-Wise function has shown itself to 
be highly amenable to lattice calculations. Most 
of the work to date has involved extracting esti- 
mates from semi-leptonic decay matrix elements 
in quenched lattice QCD at around the charm 
quark mass. Statistical errors appear to be un- 
der control and comparable in size to experimen- 
tal errors. The slope at zero recoil, — p 2 , is not 
yet precisely determined, but seems to be around 
— 1 .2 ± 50%. This uncertainty does not greatly af- 
fect the extraction of V c j and a value of 0.04±20% 
is obtained by using the lattice measurements to 
extrapolate experimental data. 

For the form factors h + and Ha 1 , which are pro- 
tected by Luke's Theorem, there are indications 



that heavy-quark symmetry holds at heavy-quark 
masses around charm and light-quark masses 
around strange. In contrast, large 0(1 /mg) cor- 
rections are observed in the form factors Zi_ and 
hy , which are not so protected. Lattice HQET 
offers a good handle on these power corrections 
by establishing the tuq — > oo limit directly. 

There is scope for substantial improvements to 
the analyses reported this year, and the possibil- 
ity of extending the calculations to semi-leptonic 
decays of the Aj. Lattice calculations of the Isgur- 
Wise function are off to an exciting start. 
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